Bound states of magnons in the S=l/2 quantum spin ladder. 
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The excitation spectrum of the two-chain Heisenberg ladder with antiferromagnetic interactions 
is studied. Our approach is based on the description of the excitations as triplets above a strong- 
coupling singlet ground state. The quasiparticle spectrum is calculated by treating the excitations 
as a dilute Bose gas with infinite on-site repulsion. Additional singlet (S = 0) and triplet (S = 1) 
excitations are found as quasiparticle bound states. Thus it is demonstrated that the spectrum 
consists of one elementary triplet, one composite triplet, and one composite singlet. 
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Quantum spin ladders have recently become a subject 
of considerable interest, mainly because of their relevance 
to a variety of quasi one-dimensional materials Q. The 
ground state properties and the quasiparticle spectrum 
of the two- leg S — 1/2 spin ladder have been studied 
by a number of techniques, including bosonization [^J, 
exact diagonalizations HL density matrix renormaliza- 
tion group calculations [H and strong-coupling expan- 
sions Js] 0| . The picture emerging from these studies is 
that of a gapped excitation spectrum and thus disordered 
ground state for any value of the inter-chain coupling Jj_ . 
Quite recently a state with energy above the one-particle 
excitation but below the two-particle continuum was ob- 
served in (VO)2P2C>7, believed to have the geometry of 
a ladder || . Subsequently this material was found to be 
better described as a dimerized chain || which in turn 
triggered theoretical work on two-particle bound states 
in this model [jioj] . The above work has raised the impor- 
tant issue of possible low-lying two-particle excitations in 
quasi one-dimensional models with a gapped spectrum. 

In the present Letter we develop a formalism to de- 
scribe two-particle bound states in dimerized quantum 
spin systems and consider the spin ladder as an example. 
We find that a scalar (S=0) and vector (S=l) bound 
state exist below the continuum and we calculate their 
properties. 

Consider the Hamiltonian of two coupled S = 1/2 
chains (spin ladder): 



H 



E 



{ JSi.Si+i + JS{.S< +1 + JA.S',} , (1) 



where the intra-chain (J) and the inter-chain (J_l) in- 
teractions are assumed antiferromagnetic J, J±_ > 0. At 
large J± 3> J the ground state consists of inter-chain 
spin singlets \GS) = |1,0)|2,0)|3,0)..., where \i,0) = 
775 [I T)i| - I Oil T);]- Since each singlet can be ex- 
cited into a triplet state it is natural to introduce a cre- 
ation operator t ai for this excitation: 



\h a ) =*Ll*>°)i a = x,y,z. 



(2) 



The representation of the spin operators in terms of t ai 
was introduced by Sachdev and Bhatt [pd|, and after its 
application to (1) one finds 



z,a,/3 

+ 2 (*ai*^i+l*/9**«*+l - 4*ai+l*0**0*+l) j • ( 3 ) 

In addition, the hard core constraint it,- it,, = has to 
be enforced on every site |11[. This exclusion of double 
occupancy reflects the quantization of spin and ensures 
the uniqueness of the mapping from (1) to (3). At the 
quadratic level the Hamiltonian (||) can be diagonalized 
by a combination of Fourier and Bogoliubov transforma- 

This gives the excitation 
J± + J cos k and 



VktL 



A%- Bl, where A 



tions t a k = v-ktak 
spectrum: u\ 
£>fc = J cos k. 

In order to address the problem of bound states of two 
triplets a reliable description of the one-particle spectrum 
is needed. We find, in agreement with previous work 
Plpl, that the effect of the quartic terms in (3) on the 
triplet spectrum is small and therefore we proceed by 
treating these terms in mean field theory. This is equiva- 
lent to taking into account only one-loop diagrams (first 
order in J). These diagrams lead to the renormalization: 

A k = Jl + J(l + 2/i) cos k, B k = J(l - 2gi) cos k, (4) 



where 



h =< *o»*crf+l >= N 1 J2q V q COS 1 and 
.91 =< t ai t m >= N^ 1 J2 q u q v q cos 1- 

The dominant contribution to the spectrum renormal- 
ization is related to the hard core condition. Previous 
treatments have used mean-field approximations [|6|JTl|] 
and are essentially uncontrolled, especially for a quasi- 
1D system. To deal with the constraint we will use the 
diagrammatic approach developed by us in Rcf. @. An 
infinite on-site repulsion is introduced in this approach 
in order to forbid double occupancy: 



Hi 



(5) 



Since the interaction is infinite, the exact scattering am- 
plitude T a f3 n $(K) = T(K)5 a ~ ( 6f3s, K = (k,u>), for the 
triplets has to be found. This quantity satisfies the 
Bethe-Salpeter equation, shown in Fig. 1(a) and depends 



1 



on the total energy and momentum of the incoming par- 
ticles K = K\ + K2- The interaction (5) is local and 
non-retarded which allows for an analytic solution of the 
equation for T: 



d 2 Q 
(2tt)2 

-T 



G{Q)G{K-Q) 



N 



LU — LU a 



(6) 



Here G(Q) is the normal Green's function (GF) G(k, t) = 
—i < T{tk a {t)t\ a (0)) > and the Bogoliubov coefficients 
u k' v k = ±1/2 + Ak/2u>k- The basic approximation made 
in the derivation of T(K) is the neglect of all anoma- 
lous scattering vertices, which are present in the theory 
due to the existence of anomalous GF's, G a (k, t) = —i < 
T{t ] _ ka {t)t\ a (0)) >. Our crucial observation (l| is that 
all anomalous contributions are suppressed by a small 
parameter which is present in the theory - the density 
of triplet excitations n = J2 a < *L*£« > ~ SN^ 1 J2 q v q- 
We find that n w 0.1 (Jl/J = 2), n « 0.25 (Jj_/J = 1) 
and it generally increases as J± decreases. Since summa- 
tion of ladders with anomalous GF's brings additional 
powers of v q into T, their contribution is small compared 
to the dominant one of Eq.(6). For consistency we also 
neglect the second term in Eq.(6). Thus the triplet ex- 
citations can be viewed as a strongly-interacting dilute 
Bose gas and Eq.(6) as the first term in an expansion in 
powers of the gas parameter n 12 1. 

The self-energy, corresponding to the scattering ampli- 
tude T is given by the diagrams in Fig. 1(b): 



nk^) = -Y,v 2 q T{k + q^-u q ). (7) 



In order to find the renormalized spectrum, one has to 
solve the coupled Dyson's equations for the normal and 
anomalous GF's. After separating the result for the nor- 
mal GF into a quasiparticle contribution and incoherent 
background, we find [ujfl: 



G{k,u) 



ZkUl 



ZkV? 



— ilk + iS lu + fife — iS 



Gi 



(8) 



The renormalized triplet spectrum and the renormaliza- 
tion constant are: 



Q k = Z k ^[A k + ^(kM 2 -B 2 k , 
8lu 



(9) 



V = 1 



j=0 



The renormalized Bogoliubov coefficients in (0) are: 



Uk ' Vk 2 + 2~n k ' 



Equations (J6|,[7jj0j|10|) have to be solved self-consistently 
for S(fc,0) and Z k . From Eq.(||) it also follows that one 
has to replace u k -> V^C/ fc , v k ^/Z^V k in (4), (j|) 
and (Q) (see also Eq.(14) below). 

Let us demonstrate how this approach works in the 
strong-coupling limit J± ^> J. To first order in J/J±, 
A k = Jj_ + J cos k and B k = J cos k. This leads to uj k w 
A k = J± + Jcosfc, u k w 1, Vk w — (J/2J±) cosk and 
h = 0, 9i = -J/4J_l. Substitution into (4), (g) and (g) 



gives 



r(fe,w)= 2J ± - lu, 

T,{k,uj)= {J/J ± ) 2 (3J ± -lj)/2. 



(11) 



Then from Eq.(|) we find the quasiparticle residue Z 



1 — (1/2)(J/J_i_) 2 and the dispersion 

n k j ± 3J j 

— = h cos k H cos 2k. 

J J 4Ji 4J, 



(12) 



The result (12) agrees with that obtained by direct 
J/J± expansion [ fl3[ to this order. For arbitrary J± a 
self-consistent numerical solution of Eqs.(||J^,g,|irj) is re- 
quired. The triplet excitation spectrum obtained from 
this solution for J±/J = 2 is shown in Fig. 3. For com- 
parison the dispersion obtained by 8-th order dimer series 
expansion [g] is also plotted. The agreement between the 
two curves is excellent which reflects the smallness of the 
triplet density n w 0.1. However for smaller values of 
the inter-chain coupling, e.g. J± = J the disagreement 
between our and the numerical results is as large as 20% 
at k — which reflects a transition into a regime, domi- 
nated by very strong quantum fluctuations. Since a more 
refined analysis is required in this regime from now on we 
will concentrate on the case J±/J — 2. 

The quartic interaction in the Hamiltonian (|^) leads to 
attraction between two triplet excitations. We will show 
that the attraction is strong enough to form a singlet 
(S=0) and a triplet (S=l) bound state. Let us consider 
the scattering of two triplets: q\a + (72/? — * 937 + <?4<5 and 
introduce the total (Q) and relative (q) momentum of 
the pair q x = Q/2 + q, q 2 = Q/2 - q, q 3 = Q/2 + p, and 
<Z4 = 0/2 — V- The bare (Born) scattering amplitude is 
(see Fig.2(a)): 



M af3n s = J {SayfipS - SapSys) COs(<? + p) + 
J (SasSpy - SapSys) COs(q - p) + 



(13) 



The J and the U terms arise from the quartic interaction 
in (3) and the constraint (||) respectively. We also have to 
take into account that the triplet excitation differs from 
the bare one due to the Bogoliubov transformation and 
the quasiparticle residue. Therefore we have to make the 
substitution: 



(10) M a p^s 



(14) 
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The bound state satisfies the Bethe-Salpeter equation 
for the poles of the exact scattering amplitude M . This 
equation is presented graphically in Fig. 2(b) and has the 
form El: 



[E Q - fi Q/2+q - n Q /2- g ] i>{q) 



dp 
2^ 



M(Q,q,p)iP(p). 

(15) 



Here M(Q,q,p) is the scattering amplitude in the ap- 
propriate channel (see Eqs. (16,19) below), Eq is the en- 
ergy of the bound state and ip(q) is the two-particle 
wave function. Let us introduce the minimum en- 
ergy for two excitations with given total momentum 
(lower edge of the two-particle continuum) Eq = 

min g {^Q/2+q + ^Q/2-g}- 11 a bound state exists then 
its energy is lower than the continuum Eq < Eq. The 
binding energy is defined as 6q = Eq — Eq > 0. 

In the singlet (S=0) channel the scattering amplitude 

is: 

M (0) = -SapS^sMafas = -4Jcosqcosp + 2U. (16) 
o 

First we consider the strong-coupling limit J± J. In 
this limit Z q — U q — 1 and fl q = J± + Jcosq. The 
edge of the continuum is Eq — 2J± — 2J\ cosQ/2| and 
therefore the Bethe-Salpeter equation (|||) is reduced to 
the simple form 



2JC Q (l + cosg) i/>(q) 



2 J cosq J ^1 cosptp(p) - U J 



ip(p). 



(17) 



We have introduced here the notation Cq — |cosQ/2|. 
The solution of (17) is: 



^ \q) = j2(l-C%) 



£ (0) 



cos q + Cq 



i0) 



/J + 2C* Q (1 -cosg) 



(18) 



We stress that the infinite on-site repulsion enforces the 
condition J dpip(p) — which means that the bound 
state is d-wave like. Thus we see that in the strong- 
coupling limit a singlet bound state always exists. At 
J± = 2 J Eq.(15) with the substitutions (16,14) has to be 
solved numerically and the result is presented in Fig. 3. 
We find that for k < 2n/5 the binding energy is practi- 
cally zero in this case. 

In the triplet (S=l) channel the scattering amplitude 

is: 



M* 1 



,7<5 



-2 J sing ship. 



(19) 



In this formula there is no summation over the index /1 
which gives the spin of the bound state. By solving the 



Bethe-Salpctcr equation in the limit J± J we obtain 
for the wave-function and the binding energy: 



smq 



(20) 



V (1) (<z) = x/i/2-2C2 — 

V eg7 J + 2C Q (l+cos<7) 

e Q 1] =2J(l/2-C Q ) 2 , C Q < 1/2. 

For Cq > 1/2 we find = which means that the 
triplet bound state only exists for momenta k > Q c — 
2ir/3 (in the strong-coupling limit). At J± = 2 J the 
numerical solution of Eq.(15), plotted in Fig. 3. shows 
that the bound state exists down to k w 7r/2. We find 
that the coupling in the triplet channel is weaker than 
the one in the singlet channel. 

The size of the bound state is determined by the spatial 
extent of the bound state wave function: 



R r 



< r z > 



(diP(q)\ 2 dq 
V dq J 2tt 



1/2 



(21) 



In the strong-coupling limit we find, by substituting 
Eqs. (18,20): 

_((i+c^(i-c Q r\ s = o 

^ s -\(l + 4C§)V2 (1 _ 4a 2 ) -i ) s = 1 W 

in units of the lattice spacing. As expected the size of 
the bound state increases with decreasing binding energy 



and near the threshold R r 



(e) 



-1/2 



e — > 0, as can be 



seen from Eqs. (18, 20, 22). The self-consistent solution for 
J± = 2J is shown on Fig. 4. Except for momenta, very 
close to the threshold, the bound states in both channels 
have typical sizes of a few lattice spacings. 

Finally, there is no bound state in the tensor (S = 2) 
channel. This can be seen from the expression for the 
scattering amplitude in this case MP) = 2J cosq cos p + 
2U which corresponds to repulsion. 

In conclusion, we have studied two-particle bound 
states in the two-leg antiferromagnetic Heisenberg lad- 
der and have shown that a singlet bound state is always 
present, while a triplet one exists only in a limited range 
of momenta. In order to find the one-particle spectrum 
we have used a new diagrammatic method to take into 
account the hard-core constraint by treating the excita- 
tions as a dilute Bose gas. Our approach is very general 
and applicable to any spin model for which the excita- 
tions can be described as triplets above a strong-coupling 
singlet ground state. We can claim that bound states 
are present practically in any system with dimcrization. 
By using the techniques described in this Letter we have 
also proven the existence of bound states in the two-layer 
Heisenberg model. Our general picture for the structure 
of the excitation spectrum fits very well into the recently 
suggested strong analogy between quantum chromody- 
namics and quantum antiferromagnetism Eq |. One can 
consider the elementary on-site spin 1/2 as a quark and 



3 



the set of triplet and singlet excitations as vector and 
scalar mesons. From this analogy one can expect that 
there are even more complex excitations, such as many- 
particle bound states, in quantum spin models of this 
type @. 
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FIG. 3. The excitation spectrum of the ladder for J± = 2 J. 
The solid lines are the elementary triplet spectrum (lower 
curve at k = n) and the lower edge of the two-triplet con- 
tinuum (upper curve). The dots are numerical results ob- 
tained by dimer series expansions [5]. The dashed and the 
dot-dashed lines represent respectively the triplet and singlet 
bound states. 

FIG. 4. The size R rm s of the bound states. The vertical 
dashed lines represent the points where the binding energy 
vanishes. 
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FIG. 1. (a) Equation for the scattering amplitude T. (b) 
Diagrams for the self-energy, corresponding to V. 



FIG. 2. (a) The bare (Born) scattering amplitude M. (b) 
The Bethe-Salpeter equation for the poles of the exact scat- 
tering amplitude M. 
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